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RESUMEN 

Este documento describe .. 

3 ABSTRACT 

04 This paper summarizes the known physics of turbulent jets observed in labora- 
^ tory experiments. The formula, which gives the power released in turbulence 

^4 describes the concentration of turbulence/relativistic particles in each point of 

the astrophysical jets. The same expression is also used to analyze the power 
released in turbulence in the case of pipe and non Newtonian fluids. Through 
an integral operation it is possible to deduce the intensity of synchrotron ra- 
diation for a profile perpendicular or not to a straight jet , a 2D map for a 
perpendicular , randomly oriented straight jet as well as a 2D map of complex 
trajectories such as NCC4061 and 3C31. Presented here is a simulation of the 
spectral index in brightness of 3C273 as well as a 2D map of the degree of 
linear polarization. The Sobcl operator is applied to the theoretical 2D maps 
^ of straight perpendicular jets. 
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^ 1. INTRODUCTION 

> 

The physical mechanism that leads to the formation of the image in extra- 
ct galactic radio-sources can be the same which explains the physics of such 

objects ; this is called the unique model. The unique model can be split in 

two 

m 

1. The extremely relativistic, mono-energetic e±-pairs of high bulk Lorcntz 
^ factor, see Kundt (1979); Kundt & Copal Krishna (2004); Kundt (2006). 

• • 

^ 2. On modeling the jets as intrinsically symmetrical, relativistic, decelerat- 

*k> ing flows it is possible to simulate the radio-map of 3C31, 0326-1-39,1553-1-24 

, NGC 315 and 3C296, sec Laing & Bridle (2002a); Canvin & Laing 
<^ (2004); Laing et al. (2005, 2006); Worrall et al. (2007). 

The transport processes are an open rather than a well established field of 
research. Following is a brief review of some approaches. 

1. The Kelvin-Helmholtz instabilities trigger a fluid turbulent cascade. The 
electrons are accelerated from the low wavelength of the turbulent spec- 
trum , see Ferrari et al. (1979); Eilek et al. (2003). 
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2. The laws of blob containment , see Zaninetti (1989) , or the 3D ran- 
dom walk , sec Zaninetti (1999) , produce simulated maps of radiation 
intensity from extra-galactic radio-sources. 

3. Plasma instabilities (e.g., Buneman, Weibel, and other two-stream in- 
stabilities) created in collision-less shocks may be responsible for particle 
(electron, positron) acceleration, see Nishikawa et al. (2006) . 

The canonical approach to the turbulence is through the Navier-Stokes 
equations 

^ = lvp + .V^U . (I) 

where U represents the velocity vector , is the fluid particle acceleration 
, v is the kinematic viscosity, p the density and p the pressure. By applying 
the Reynolds decomposition to the previous equations 

U(x,t) = (U(x,t)) + u(x,t) , (2) 

where u(x,t) represents the fluctuation, we obtain the equation of the mean 
momentum or Reynolds equations 

^ = -.V^(t/,)-^-^f^ , (3) 

Ot OXi p OXi 

where |^ represents the mean substantial derivative. The difference between 
the Reynolds and the Navier-Stokes equations is represented by the Reynolds 
stresses , ^^^^"^^ . The following reviews three different models which adopt 
the turbulence in astrophysical jets. 

1. A complex approach to the turbulence is through the ensemble by Favre 
(1969) , see Bicknell (1984) , where the density p is expressed as p + p' 
and the velocity Vi as Vi+v'^ where {i = 1, 2, 3). From this double decom- 
position it is possible to obtain four equations for mass , momentum , 
turbulent kinetic energy and heat. The behavior of the magnetic fleld is 
obtained by taking into consideration the ensemble average of the MHD 
equations for the magnetic field Bi in the infinite conductivity limit. 
Once a law for the pressure of the atmosphere is considered, it is possi- 
ble to obtain a law for the velocity , the electron distribution function 
, the magnetic field dependence and the density. The above equations 
allow us to deduce a law for the brightness along the jet I'^ once an adi- 
mensional opening angle , a' = (p/ipo , is introduced. The brightness law 
is of the type II oc with b [0.35 — 0.78] . In this model the opening 
angle is variable , and the application was done to 3C31 which shows 
a variable opening angle starting from 12.3° , reaching a maximum of 
16.6° in the middle and re-coUimating in the outer regions at 12.6° , see 
Table 1 in Bridle et al. (1980). Another application of this model was 
the study of the physical state of the interstellar medium in the optical 
counterpart NGC 383. 
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2. A compressible k-e turbulence model was used by Falle (1994) to study 

the effect of turbulence on the propagation of fluid jets. In this model a 
connection between the opening angle of the jet 6 and the Mach number 
, M , is 

M=\ . (4) 

Through a numerical simulation the z-component of velocity along a cut 
perpendicular to the z-axis has been evaluated. 

3. The near infrared images of HH 110 jet (Herbig Haro) were interpreted as 

due to low velocity shocks produced by turbulent processes, see Noriega- 
Crespo et al. (1996). The spatial intensity distribution of H2 , Ha and 
[•S*//] 617/31 perpendicular to the flow axis and along the cross section of 
knots in HHllO has a behavior that can be approximated by a Gaussian 
distribution , see Figure 4 in Noriega- Crespo et al. (1996). In one case , 
Ha in knot P of HHllO, it is possible to see a bump near the maximum 
of the intensity in the transversal direction. This can be considered the 
first observational evidence of a physical effect which will , later on , 
be called " valley on the top" . The results were explained by Noriega- 
Crespo et al. (1996) adopting a turbulent mixing layer that followed a 
Couette flow. 

The previous approaches do not solve the following questions 

• Can the physics of turbulent jets , which are observed in the laboratory, 
be applied to extra-galactic radio-sources ? 

• Is it possible to build the image of a straight turbulent jet which is emit- 
ting synchrotron radiation from the law of dissipation of turbulence ? 

• Could the observed cuts in non-thermal radiative flux, parametrized as 
a function of the distance , be compared with theoretical cuts ? 

• Could the 2D maps in brightness that characterize the complex mor- 
phologies associated with radio-galaxies be simulated ? 

• Is it possible to simulate the variation of the spectral index in brightness 
of the non-thermal emission ? 

The previous cited papers on turbulence solve some of the problems con- 
nected with the application of turbulence to astrophysical jets but leave other 
problems open and Table 1 reports the status. 

This paper briefly introduces in Section 2 an approximate law for the con- 
version of energy of the bulk flow into the turbulent cascade and the radiative 
transfer equation. Section 3 computes the power released in the turbulent cas- 
cade for three type of fluids : turbulent fluid , pipe fluid and non- Newtonian 
fluid. Section 4 computes through an integral, the intensity profiles of the jet 
when it is directed in the direction perpendicular to the motion for the three 
types of fluids previously analyzed. 
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TABLE 1 

SYNOPTIC TABLE OF THE SOLVED AND UNSOLVED PROBLEMS IN 
THREE DIFFERENT TURBULENCE PAPERS AND HERE 



Problem 


Bicknell 1984 


Falle 1994 


Noriega — Crespo et al. 1996 


this paper 


brightness along the jet 


yes 


not 


not 


yes 


perpendicular brightness 


not 


not 


not 


yes 


image as integral operation 


not 


not 


not 


yes 


law for magnetic field 


yes 


not 


not 


not 


valley on the top 


not 


not 


yes 


yes 



Section 5 deals with an algorithm used to build the synchrotron image in 
two simple cases : the straight jet is oriented in the direction perpendiculax 
to the observer or has a random direction. 

Section 6 presents the images of the radio-jets which present bending and 
wiggling , such as NGC4061 and 3C31. This Section also contains a simulation 
of the spectral index in brightness of turbulent astrophysical jets. 

2. PRELIMINARIES TO THE RADIO-IMAGES 

Section 2.1 analyzes a possible mechanism for conversion of energy from 
the flux of kinetic energy down to the turbulence. The matrix equation that 
implements the radiative transfer equation is introduced in Section 2.2 . 

2.1. The efficiency of energy conversion 

The total power , Q , released in a turbulent cascade of the Kolmogorov 
type is, see Pelletier & Zaninetti (1984) 

Q ~ IkhPSt , (5) 

where jkh = lad^ is the growth rate of K-H instabilities , p is the matter 

density and st is the sound velocity, see Section 9.1 in Zaninetti (2007). 

The total maximum luminosity , Lt , which can be obtained for the jet in 
a given region of radius aj and length aj is 

Lt = -Ka^ajQ = 'yadPSr'^a] . (6) 

The mechanical luminosity of the jet is 



(7) 
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and therefore the efficiency of the conversion , xt , of the total available 
energy in turbulence is 

^* 2 11 

where M is the Mach number. Prom this equation it is clear that the fraction 

of the total available energy released firstly in the turbulence and after in non 
thermal particles is a small fraction of the bulk flow energy. The assumption 
made in Section 3 and Section 4 by Zaninetti (2007) in which the bulk flow 
motion was treated independently from non-thermal emission is now justified. 

2.2. Radiative transfer equation 
The transfer equation in the presence of emission only , see for exam- 
ple Rybicki & Lightman (1985), is 

f--' ■ 

where 1^ is the specific intensity , s is the line of sight , jv the emission 
coeSicient, and the index v denotes the interested frequency of emission. The 
solution to equation (9) is 

h{s) = I^so) + f 3^{si)dsi , (10) 

J So 

where /y(so) represents the intensity at s = sq , in this case zero. In the case 
of synchrotron emission, the maximum intensity of mono-energetic electrons 
in the presence of a magnetic field B in gauss is 

his) = 1.710-2^5^ f c^^^^s /""^ ^ , (11) 
J cm-' s sr Hz 

where C(s) is the number of electrons in a unit volume along the line of 

sight and B± is the component of B perpendicular to the velocity vector , 
see for example equation (5.550) in Ginsburg (1975). In the case of spatial 
dependence on it is assumed 

> = QC(s) , (12) 

where 

C,- = 1.710-2^5^ • (13) 

s sr Hz 

The intensity in the non-homogeneous case is 

/^(s) = Cj I C{si)dsi . (14) 

J sa 

The increase in brightness is proportional to the concentration integrated 
along the line of sight. In the Monte Carlo experiments here analyzed, the 
concentration is memorized on the three-dimensional grid C and the intensity 
is 

I{i,j)=Y,5xC{i,3,k) , (15) 
fc 

where 5 is the spatial interval between the values of concentration and the 
sum is performed over the interval of existence of index k. 
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3. THREE TYPES OF FLUIDS 

The cases of shear layer produced by a turbulent jet, by a pipe and by a 
non-Newtonian fluid are now analyzed. 

3.1. The turbulent jet 

The theory of turbulent round jets can be found in different textbooks. 
The more important formulas are now reviewed as extracted from chapter V 
by Pope (2000) ; similar results can be found in Bird et al. (2002) and 
in Schlichting et al. (2004) . We start with the centerline velocity Uo{z) , 
equation (5.6) in Pope (2000) , as measured in the laboratory experiments : 

^ = A (16) 

C/i z/d ' ^ ' 

here z denotes the main direction , d is the diameter of the nozzle, B is 
a constant derived in the laboratory that takes the value 5.8, and Ui is the 
initial jet velocity. The solution of the mean velocity <U > , equation (5.100) 
in Pope (2000) , along the main direction is 

SttTl^T 1 

<[/>=f/i^^— , 17 

z (1 + aT?? ) 

where = ^ , r is the radius of the jet at z, ax is a constant and vt is the 
turbulent viscosity. The viscosity , equation (5.104) in Pope (2000) , is 

UT = ^ , (18) 

and ax, equation (5.18) in Pope (2000), is 

ar=^-^ . (19) 

where S will be later defined. The production of turbulent kinetic energy in 
the boundary layer approximation , equation (5.145) in Pope (2000) , is 

where y is a Cartesian coordinate that can be identified with r. The flow rate 
of mass m(z) is , see equation (5.68) in Pope (2000) , 

/■oo 

miz) = 2wph/2{h/2iz)Uo{z))x / ^/(Orf^ , (21) 

Jo 



where 



and 



1 

(1 + Aey 



m = jr^^ , (23) 
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TABLE 2 

PARAMETERS OF THE TURBULENT JET WHEN k= 0.54 . 



a[rad\ 


a[degree] 


S 


B 




0.087 


5 


0.043 


12.37 


0.013 


0.185 


10.64 


0.093 


5.79 


0.028 


0.343 


20 


0.17 


3.06 


0.053 


0.523 


30 


0.26 


2.01 


0.08 



where A is a constant that will be later defined and 61/2 is the value of the 

radius at which the velocity is half of the centerhne value. The jet draws 
matter from the surrounding mass of fluid. Hence, the mass of fluid carried by 
the jet increases with the distance from the source. The previous formulas are 
exactly the same as in Pope (2000) ; we now continue toward the astrophywical 
applications. The quantity S is connected with the opening angle a through 
the following relationship 

5 = tan I . (24) 

The self-similar solution for the velocity , equation (17) , can be re-expressed 
introducing the half width z = bi/2/S 

where A = \/2 — 1 . From the previous formula the universal scaling of the 
profile in velocity is now clear. From a careful inspection of the previous 
formula it is clear that the variable z should be expressed in d units (the 
nozzle's diameter) in order to reproduce the laboratory results. In doing so 
we should find the constant k that allows us to deduce B 

B = kx SoTfT ■ (26) 

Table 2 reports a set of S , S and vt for different opening angles a . The 
assumption here used is that k is the same for different angles. The velocity 
expressed in these practical units is 

BUi 1 



<U>=^ ,, . r ^2^2 " (27) 



This formula can be used for z expressed in d-units when z > B and 
Figure 1 reports the field of velocity for a laboratory jet. 

The first derivative of the profile in velocity as given by formula (27) with 
respect to the radius is 

^ <U>=UiX ^- '-^ = . 28 

tan (1/2 a) (^2 - l) ^ (6,// + A r^y 



8 



ZANINETTI 



minimum= 0,0004 maximum= 0.2998 




10 20 30 40 



z 

Fig. 1. Theoretical 2D map of the velocity ,<(/>, given in units of the sound 
speed when the initial Mach number is 0.3. The length considered is z/d — 44 , k 
=0.54, A = 0.414 and adeg = 10.64 . 
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The production of turbulent kinetic energy is 



2 ^/2-1^4fc26 /2^V2 

V = vtUI X ^- ^ '-^ ^ , (29) 

tan (1/2 a) {^2 ~ l) z2 (^^^^2 ^ ^ ^2)'' ^^^^4 

and z should be expressed in ti-units. It is interesting to note that the maxi- 
mum of , is at 

r= ■ (30) 

A quantity that is measured in the laboratory is the turbulence intensity 

Tu=^ , (31) 



Tu 



where Urms represents the standard deviation of the velocity measured in the 
laboratory and < U > the averaged velocity. This quantity is w 25% at 
r = and rises with increasing r , see Pope (2000) and Hussein et al. (1994). 
Recent results obtained using the PIV (Particle Image Velocimetry) are able 
to produce sophisticated measurements , see Figure 2 , extracted from Bravo 
(2006) , who reports a 2D map of Tu . A possible theoretical explanation for 
Tu can be found from the only combination for Urms derived from V which 
produces an adimensional number when divided hy < U > 



'vT d<U> 




Fig. 3. Theoretical 2D map of the turbulence intensity , Tu , when Tumax ~ 30%. 
The length considered is z/d = 44 , =0.54, A = 0.414 p =1 , and ad^g = 10.64 . 



where p is the density. Using the practical units 



Tu= V^K^v^-+^-) y_ ' (33) 



Figure 3 reports the theoretical 2D map of Tu as given by equation (33) 
normalized to a maximum value Tumax ^ given by laboratory experiments. 
When the flow rate of mass is expressed in these practical units equation (21) 
becomes 

^. ^ ^pz(tan(l/2a))^t/,rf 

In this formula both z and d are present and it is possible to speak of physical 
units, but z > d . 



3.2. Pipe fluids 



In smooth circular tubes the averaged velocity for turbulent flow in the 
direction perpendicular to the motion vj_ is often parameterized through the 
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Fig. 4. Mean velocity profile ( in relative units ) vs. jet radius (in relative units) 
when n=6 ( full line ), n=7 ( dashed ), n=8 ( dot-dash-dot-dash ), n=9 ( dotted ) 
and n=10( dash-dot-dot-dot ). 



following equation 

^ = (l-^)i, (35) 

where r is the distance from the jet central axis , a is the jet radius, v±^max is 
the maximum radial velocity and n is an integer, see for example Bird et al. 
(2002). Figure 4 reports the mean velocity profile for different values of n. 
On introducing AR = a — r , the following is found: 

^ = (^)4^ (36) 

Further on n is connected to the Reynolds number ,Re, through the following 
relationship 

n = 7V/iVr(10"«'= X R'"'--), (37) 

where NINT represents the nearest integer. The couple (aiie,^fle) is easily 
found once the experimental correspondence n — > Re is provided , see for 
example Table 3 or Bird et al. (2002). On applying the least square fitting 
procedure , a^e =0.561 and hf>e =0.0582. Once a value of — is chosen 
we can plot ^ as a function of the Reynolds number, see Figure 5 . Table 3 
summarizes the equivalence between n, Re and when — is fixed . The 
first derivative of the profile of perpendicular velocity as given by formula (35) 
is then computed 
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a: 
< 

3 - 




3 - 

ftclOf" 10" IJbcIO™ 2x10'" iOxlO" 3x10*" ISxIO' 

Re 

Fig. 5. ^ vs. Reynolds number when — — =0.9 . 



TABLE 3 

EQUIVALENCE BETWEEN PARAMETERS WHEN ^.^^^^ =0.9 



Re ^ 



6 4 X 10^ 0.53 

7 1.1 X 10^ 0.47 

8 1.1 X IQS 0.43 

9 3.2 X 10*^ 0.38 



JET IMAGE 



13 



By analogy with formula (20) a power released in the cascade as 

rn = u„i^f , (39) 
is introduced where the index n stands for pipe fluid. We therefore have 

(- 1 \ 2 

rn= \ V ^ • (40) 

(—a + r) 

3.3. Non-Newtonian fluid 
In a non- Newtonian fluid the shear stress , Tyr , is often parameterized as 

dr' 

where kn is a constant , u the velocity , r the distance from the center and m 
an integer. The profile in velocity results to be 

i*i = i-(r)=s^ , (42) 

where a is the radius of the pipe; see Problem 7-1 in The Staff of REA (1985) 
or Section 12.3 in Hughes & Brighton (1967). Also here we compute the first 
derivative of the profile of perpendicular velocity 

d u(r) _ {l)"^{m+l) 



dr Umax mr 
The power released in the cascade is 

,du{r) 2 



(43) 



T^non — ^noni ^ ) ; (44) 

or 

where the index non stands for non-Newtonian fluid. We therefore have 

fnon iumax{'i)~\ {m + if 
Vnon = ^ Yir^ • (45) 

4. PROFILES IN THE PERPENDICULAR DIRECTION 
The following explores how the profile in intensity of a jet intersected by 

a plane perpendicular to the direction of motion should appear. 

Therefore, a classic Sobel filter is applied in order to find gradients in the 

profiles. 

It is assumed that the synchrotron-emissivity scales as the power released 
in turbulent kinetic energy, see equation (29), 

e^V . (46) 

The cross section of the jet is reported in Figure 6. Due to the additivity 
of the synchrotron emission along the line of sight, an integral operation is 
performed in order to obtain the intensity of emission 

I{y) = 2 X e{r)dx , (47) 

Jo 

with r = \/ -\- y'^ and a representing the jet radius. 
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Fig. 6. The conical source is represented through a section perpendicular to the jet 
axis. The observer is situated along the x direction, one line of sight is indicated. 




Fig. 7. Intensity of radiation I{y) across the turbulent jet when z/d = 44 , k =0.54, 
A = 0.414 and aaeg = 10-64 . 

4.1. Profiles of the turbulent jet 
The intensity of emission according to formula (29) is 

Jo {bi/s'^ + Ax'^ + A y^y (\/2 - l) tan (1/2 a) 

(48) 

This integral has an analytical solution but it is complicated and therefore 
Figure 7 only shows the numerical integration. The numerical integration 
gives a characteristic shape on the top of the blob which is called the " valley on 
the top" and the observational counterpart of this new effect will be discussed 
in Section 4.4. The maximum of this integral is at the point y ~ 0.491356i/2 
and the value of intensity at the maximum is 1.18 times the value at the point 
2/ = ( the center of the jet). 

4.2. Profiles of the pipe fluids 
The intensity of emission in pipe fluids according to formula (40) is 



i{y) ~ / 

Jo 



max 



dx . (49) 



rfi [a — ^Jx"^ + 

The previous integral has an analytical solution in the case n = 1 

J(y) = 2 ''"^^-'"°"y"'~^' , (50) 
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Fig. 8. Intensity of radiation I{y) across the pipe fluid jet, when = 1 , a = 1 , 
v±,max = 1 and n = 1 . 

which is reported in Figure 8. 

When n 1 the integral cannot be solved in terms of any of the standard 
mathematical functions 

and Figure 9 reports the numerical integration . 

4.3. Profiles of non-Newtonian fluid 
The intensity of emission in a non-Newtonian fluid accovding to formula (45) 

is 



i{y) 



m? {x^ + y2) 



-dx 



(51) 



Inserting Vnon=l , a = 1 and Umax=^ for the sake of simplicity, the previous 
integral has the analytical solution 

i{v) = 

m+1 , -1-1 1 — iP" « 

2(t/2) yr^2i^i(l/2,l-"^;3/2; ^)(rn+l)y-^m-^ + 

m + l ^ _l_ 1 1 9(2 

2(2/2) - 71^7 2^,(1/2, 1-"^;^^; 3/2; -^)(m+l)m-V' (52) 



m 



where 2Fi{a, b; c; z) is a regularized hyper-geometric function, see Abramowitz 
& Stegun (1965), von Seggern (1992) and Thompson (1997). Figure 10 reports 
the results of the numerical integration. 
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Fig. 10. Intensity of radiation I{y) across the non-Newtonian jet normalized to 1 

when ^'nori=l , a = 1 and Umax~i, and m = 1 (full lino ), m = 2 (dashed ), m = 3 
(dot-dash-dot-dash), m = 4 (dotted) and m = 5 (dash-dot-dot-dot) . 
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Fig. 11. X-intensity of knot D-E in M87 (empty stars ) and theoretical intensity 
of radiation for turbulent fluids I{y) (dotted line) when z/d = 44 , k =0.54, A = 
0.414 and adeg = 10.64 . 

4.4. Comparison with observations 

This paragraph starts by checking the profile in the direction perpendicular 
to the jet of M87 , in particular to the X-component of knot D-E in Figure 2 
by Perlman & Wilson (2005). The data of the observations are reported 
in Figure 11 together with the turbulent jet profile scaled in order so that 
the value of half theoretical and observed intensity are approximately equal. 
Figures 12 and 13 report the comparison of pipe-fluid and non-Newtonian 
fluid respectively. 

Another interesting situation is represented by the four blobs in 3C273 as 
imaged at the frequency of 5GHz by the VSOP (VLBI ( Very Long Baseline 
Interferometry) Space Observatory Programme) . Here the focus is on blob A 
which presents a symmetry around the center of the blob and on the "valley 
on the top" . 

The comparison between observed and theoretical-turbulent "valleys on 
the top" , see Figure 14, shows that the theoretical depression is bigger than 
the one observed. In the case of non-Newtonian fluid m = 5 produces a perfect 
agreement between observed and theoretical depressions for blob A of 3C273, 
see Figure 16. 

4.5. Sobel Filter 

The Sobel filter , see Pratt (1991) , computes the spatial gradient of the in- 
tensity of radio-galaxies , see , for example, Figure 4b and Figure 26b in Laing 
et al. (2006), which represent respectively 3C296 and 3C31. The Sobel filter 
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Fig. 12. X-intensity of knot D-E in M87 (empty stars ) and theoretical intensity of 
radiation for a pipe-fluid I{y) (dotted line) when = 1 , o = 1 , v±^rnax = 1 and 
n = l . 
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Fig. 13. X-intensity of knot D-E in M87 (empty stars ) and theoretical intensity of 
radiation for a non non-Newtonian fluid I{y) (dotted line) when i/„o„=l , a = 1 , 
Um.ax=i, and m = 5 (dotted line). 
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Fig. 14. VSOP at the frequency of 5GHz of region A in 3C273 (empty stars ) 
and theoretical intensity of radiation for turbulent fluids , I{y), (dotted line) when 
z/d = AA ,k =0.54, A = 0.414 and aaeg = 10.64 . 



is now applied in the following way . S is the the source image , Gx and 
Gy are the two images which at each point contain the horizontal and ver- 
tical first derivative of S. At each point in the image the resulting gradient 
approximations can be combined to give the gradient magnitude G 



G = a/G 



G 



y ' 



(53) 



where 



G, 



and 



10-1 

2 0-2 
10-1 

1 2 1 

-1 -2 -1 



*S 



(54) 



(55) 



where * denotes the 2-dimensional convolution operation. The gradient direc- 
tion, Q , is 

e = arctan . (56) 

Here we have implemented the Sobel filter to ID cuts in intensity of 3C273; 
Figure 15 reports the theoretical and the observed Sobel filter. 
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[mas] 



Fig. 15. Gradient of the intensity (Sobel filter) for VSOP-data at the frequency of 
5GHz of region A in 3C273 (empty stars ) and gradient of theoretical intensity (Sobel 
filter) of radiation for turbulent fluids (dotted line ) . Parameters as in Figure 14. 




Fig. 16. VSOP at the frequency of 5GHz of region A in 3C273 (empty stars ) and 

theoretical intensity of radiation for a non non- Newtonian fluid I{y) (dotted line) 
when i'non=i , a = 1 , Umax=i , and m = 5 (dotted line). 
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nninimum= —5.3405 nnaxirrium= 2.0000 
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Fig. 17. Theoretical 2D map of the decimal logarithm of the intensity of emission 
over the astrophysical length of 20kpc corresponding to the radio-galaxy NGC4061 
at 21 cm , see Figure 5 in Zaninetti (2007) . The length expressed in the nozzle's 
diameter units is z/d = 44 , k =0.54, A = 0.414 and Odeg = 10.64 . The integral 
operation is performed on a cubic grid of 400^ pixels. 

5. 2D MAPS 

An algorithm is now built , which is able to simulate the synchrotron 
emission of the jet given the basic assumption that the emissivity e scales as 
the power released in turbulence for a turbulent fluid. The integral operation 
that gives the intensity on a plane is performed through a simple addition as 
given by formula (15). The simplest case , a jet perpendicular to the line of 
sight is now analyzed. 

The symmetry here used is that the jet's main direction is perpendicular 
to the observer, see Figure 17. 

In order to clarify the interpretation of this theoretical 2D color map Fig- 
ure 18 reports the intensity profile along the main direction and the transversal 
profile as computed in the middle of the jet is shown in Figure 19. 

The previous case represents a cut along a direction that forms an angle 
of 90° with respect to the jet's main direction. When the angles are different, 
the " valley on the top" becomes asymmetric and the thickness of the cut 
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Fig. 18. Cut of the intensity of emission along the main direction at the center of 
the jet as extracted from Figure 17 . 




r [kpc] 



Fig. 19. Cut of the intensity of emission in the transversal direction , or 90° in the 
middle of the jet as extracted from Figure 17 . 
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Fig. 20. Cut of the intensity of emission in the middle of the jet at 45° (full line) 
and at 120° (dashed line) as extracted from Figure 17 . 



increases , see Figure 20 where two profiles are reported. 

The gradient magnitude of the previous image is now plotted , see Fig- 
ure 21. 



6. ASTROPHYSICAL APPLICATIONS 

The following reports the velocity , the equation of motion , the power 

released in turbulence and the flow rate of mass when the astrophysical units 
are adopted. The algorithm that allows us to build complex trajectories such 
as NGC6041 and 3C31 is then introduced. The observed spectral index typical 
of the radio-galaxies is simulated. 



6.1. Astrophysical equations 

Here shown is the astrophysical version of the previous derived equations 
for the turbulent jet. Equation (16) allows us to deduce the centerline velocity 

«o(.) = 162027 —-§1^ , (57) 

^ ^ tan (0.00872 adeg) -21 sec ^ ^ 

where adeg is the opening angle expressed in degrees, Pq is the initial velocity 
divided by the light velocity, di is the diameter of the nozzle in kpc units 

and zi is the length of the jet in kpc units. From the previous equation it is 
possible to deduce the equation of motion for a turbulent jet , 



^(i) = 57.56j iMlll kpc , (58) 

^ ' y tan (0.00872 adeg) 




Fig. 21. Sobel filter ol the intensity of the central part of NGC4061. The gradient 
magnitude is computed according to formula (53), other parameters as in Figure 17. 
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where tr = t/ (lO'') year. The radius of the turbulent jet is 

K^) = [y + ^Wtan(^)]/cpc , (59) 

where a is the opening angle expressed in radians. The power released in the 
turbulent cascade is 

e(..) = 2(v/2-l)-F5,,/-AV ^ 

tan(l/2a)(v/2-l)'(6,// + Ar2)%,// ^1 ' 

where zi is the position of the jet expressed in kpc units , A and k are defined 
in Section 3.1 . 

The flow rate of mass , see equation (34) , as expressed in these astrophys- 
ical units is 

m{z) = 8.0 10" no zi [tan (|)] ' /3o di , (61) 

where po = nom , m = 1.4 mu , no is the density expressed in particles cm~^ 
, mu is the mass of the Hydrogen and A4q is the mass of the Sun. 



6.2. Astrophysical images 

The points belonging to a trajectory , in the following TP , are now inserted 
on a cubic lattice made by pixels^ points . The intensity in each point of a 
2D grid is computed according to the following rules. 

1. In each point of a 3D grid we compute the nearest TP. 

2. The concentration on a 3D grid is then computed because we know the 
distance between the grid point and the nearest TP. 

3. The intensity of radiation is then computed according to the integral 
operation as represented by formula (15). 

A first test can be done on a radio galaxy already examined in Zaninetti 

(2007), NGC4061. The equation of motion is formula (58), derived in the 
previous paragraph. The other parameters necessary to produce bending and 
wiggling of the jet as well as the Eulerian angles which represent the point of 
view arc reported in the caption of Figure 22. 

The phenomena of the "valley on the top" can be visualized through a 
surface rendering of the 2D intensity , see Figure 23 , or through a ID cut of 
the matrix that represents the intensity , see Figure 24 . 

Another radio-galaxy that can be here simulated is 3C31 ; this radio-galaxy 
has been chosen because it does not present knots, see radio observations in 
Figure 25. This radio-galaxy has already been simulated in the first straight 
12 kpc , and a comparison between the model and the observed brightness 
distribution is reported in Figure 6 of Laing & Bridle (2002b) and in Fig- 
ure 3 of Laing et al. (2005) . The first 60 kpc are simulated in order to 
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Fig. 22. Theoretical 2D map of the decimal logarithm of the intensity of emission 
representing the continuous three-dimensional trajectory of NGC4061: the three 
Eulerian angles characterizing the point of view are $=80 ° , Q=-85 ° and ^l/^-lO °. 
The precession is characterized by the angle *l/prec= 5 ° and by the angular velocity 
^iprec= 60.00 [°/10^year]. The three Eulerian angles which characterize the jet are 
$j.= ° , ^'j= ° and 6^= 90 °. The angle of rotation of the galaxy is aG= °. 
The physical parameters characterizing the jet motion are : t7= 12.00 lOVear, di = 
0.01 kpc , /3o = 0.1 , ao= 10.64 °. The length of the jet is 20.6 kpc and the velocity 
of the galaxy Vz= 117 km/s . The integral operation is performed on cubic grid of 
500^ pixels. 



28 



ZANINETTI 




Fig. 24. Intensity of NGC4061 represented through a cut in the x-direction , pa- 
rameters as in Figure 22. 
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cover the wiggles. The 3D surface that represents the trajectory is presented 
in Figure 26, the simulated intensity distribution is reported in Figure 27 , 
Figure 28 reports the gradient magnitude computed according to the Sobel 
filter, Figure 29 reports the intensity represented through a surface rendering 
that visualizes the "valley on the top", and Figure 30 reports the intensity 
represented through a ID cut. 



6.3. Spectral index 

The collision of an electron of energy E with a cloud produces an average 
gain in energy (A£J) proportional to the second order in u/c , u is cloud 
velocity and c is light velocity, 

This process is named Fermi II after Fermi (1949). 

On introducing an average length of collision , A , the formula becomes; 

where 

It should be remembered that with a mean free path between clouds A, the 
average time between collisions with clouds , see Lang (1999) , is 

2- . (65) 
c 

The probability . P{t) , of the particle remaining in the reservoir for a 
period greater than t is, 

P{t)=e-T , (66) 

where T is the time of escape from the considered region. The hypothesis 
that energy is continuously injected in the form of rclativistic particles with 
energy Eq at the rate R, produces ( according to Burn (1975)) the following 
probability density , N , which is a function of the energy: 

with 

7=1+^ , (68) 
and T is defined in equation (64) . Equation (67) can be written as 

N{E) = K E-^ , (69) 

where K = ^J^'^+i ■ A power law spectrum in the electron energy has now 
been obtained. A power law distribution of relativistic electrons ,see Lang 
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Fig. 25. The large scale structure of 3C31 at the frequency of 1.4 GHz. The north- 
south field covers 300 kpc and the small box here modeled covers 120 kpc. This image 
was generated with data from telescopes of the National Radio Astronomy Obser- 
vatory, a National Science Foundation Facility, managed by Associated Universities. 
Inc . 
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Fig. 26. Continuous three-dimensional surface of 3C31: the three Eulerian angles 
characterizing the point of view are "1>=180 ° , 0=170 ° and \&=195 °. The precession 
is characterized by the angle ^prec= 10 ° and by the angular velocity np^ec= 54.00 
[°/10^year]. The three Eulerian angles which characterize the jet are ° , = 
° and Oj= 90 °. The angle of rotation of the galaxy is ac— °. The physical 
parameters characterizing the jet motion are : t7= 10.00 lO^year, di= 0.2 kpc , f3o 
= 0.05 , ao= 10.64 °. The length of the jet is 59.65 kpc and the velocity of the 
galaxy 11^= 292 km/s . 
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Fig. 27. Theoretical 2D map of the decimal logarithm of the intensity of emission 
representing the continuous three-dimensional trajectory of 3C31; parameters as in 
Figure 26 . The integral operation is performed on a cubic grid of 500^ pixels. 
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minimum= 0.0000 maximum= 10.0000 

I 





Fig. 28. Sobel filter of the intensity of tiie central part of 3C31. The gradient 
magnitude is computed according to formula (53), parameters as in Figure 26 . 
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Fig. 30. Intensity of 3C31 represented through a ID cut along the x-direction , 
parameters as in Figure 26 and Figure 27. 
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(1999) , corresponds to a power law frequency spectrum proportional to v"" 
, where the spectral index is 

as = --^ = -Zf^ ■ (70) 

T is now assumed to be constant and r scales with the power released in 
turbulent kinetic energy, e , as 

1 



' [<r,z)]^/b ' ^^^^ 

where 6 is a parameter that connects theory with observations. Concerning 
the observations of the spectral index along the jet we concentrate on 3C273 , 
see X-observations by Jester et al. (2006) and far-ultraviolet imaging by Jester 
et al. (2007). Table 3 in X-Ray by Jester et al. (2006) was selected as an ex- 
ample to simulate : the spectral index is 0.83 at knot Bl (internal region) 
and 1.27 at knot H2 (external region). Therefore as suggested by the theoret- 
ical formula (70) the observed spectral index increases in absolute value going 
from the inner regions to the outer regions. The theoretical spectral index 
was simulated in the following way 

1. A spectral index is generated according to formula (70) in each point of 
the 3D jet. A typical 2D cut is reported in Figure 31. 

2. The intensity at the frequency v, 1^ , is computed in each point of the 
jet according to equation (60). 

3. The intensity at a frequency one decade bigger , u x 10, is evaluated 
according to the scaling , I^xio = IvilO)"" ■ 

4. The integral operation is performed on a cubic grid of pixels^ points. 
The spectral index that comes from the two integral operations at two 
different frequencies is computed , see Figure 32. The parameter b is 
chosen in order to fit the observations. 

In this simulation , see Figure 32 , the spectral index increases in absolute 

value both in going from the internal regions to the external regions following 
the axial direction and in going from the center of the jet to the boundary along 
the perpendicular direction. This effect has been observationally mapped in 
the case of 3C296 , see for example Figure 4b in Laing et al. (2006) . 

Under the hypothesis of ultra-relativistic electrons , a power law —7 in the 
energy distribution of the electrons and a uniform magnetic field, the degree 
of linear polarization is, see equation (1.176) in Lang (1999) , 

^ (-) 

Figure 33 maps the degree of linear polarization that comes from the 2D map 
of the spectral index after a conversion of the index 7 in ofg , 




Fig. 31. Cut of the first spectral index through a plane crossing the center of the 
jet, 6=14. The length of the jet is 24fcpc , the same as 3C273. The 2D grid is made 
of 400^ pixels. 




Fig. 32. Spectral index of 3C373 as comes from the integral operation over two 
frequencies. The integral operation is performed on cubic grid of 400^ pixels. 
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. 33. Degree of linear polarization of 3C373; parameters as in Figure 32 . 
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7. SUMMARY 

Theoretical Emissivity The power released in the turbulence is com- 
puted for three type of fluids. The maximum input in turbulence is realized 
when r = 0.696i/2 , in the case that turbulent fluids are considered. 

Theoretical Intensity Profiles The profiles in intensity of synchrotron 

radiation perpendicular to the line of sight are analytically or numerically 
solved when the simplest configuration is realized : straight jets oriented in 
the direction perpendicular to the observer. The maximum of synchrotron 
radiation for turbulent fluids is at r w 0.496i/2- This means that a depression 
at the center of the jet is expected. The three types of synchrotron radiation 
here analyzed are then compared with the observed X-component of knot D-E 
in M87 and with knot A in 3C273 at the frequency of 5GHz. The confrontation 
with the data from astronomical observations should be carefully done because 
the data are often smoothed and this operation may change the structure of 
the central depression. 

Theoretical Radio maps More complex is the set up of an algorithm 
that implements the radiative transfer equation for a randomly oriented jet. 

In order to build such an image, the centerline trajectory should be made 
discrete and the distance between the point of a 3D grid and the nearest 
point of the trajectory computed. This algorithm computes the theoretical 
map in intensity of synchrotron radiation for complex morphological cases 
such as NGC4061. The depression at the center of the jet , Figure 23, can be 
considered an alternative to the cosmic double helix , see Lobanov & Zensus 
(2001); Lobanov et al. (2005). 

Theoreticcd Spectral Index The connection between power released in 
turbulence and acceleration of electrons has been attached adopting a Fermi 
mechanism with a variable velocity of acceleration. This mechanism, after the 
introduction of a parameter that connects theory and observations allows us 
to build 2D maps of the theoretical spectral index of synchrotron emission, see 
Figure 32. The simulation predicts that the spectral index should decrease ( it 
increases in absolute value) along the centerline and the transversal directions. 
As a consequence the 2D map of the degree of linear polarization is easily built. 

Theoretical Sobel Filter The new technique of the Sobel Filter can de- 
tect details in the map of radio-galaxies. The phenomena of the "valley on 
the top" can be revealed due to the gradient's change of sign in going from 
the center to the external regions in a radial direction, see Figure 15 and Fig- 
ure 21 . The presence of arcs in the Sobel-filtered images , see Figure 26b in 
Laing ct al. (2006) , can be explained by the fact that the integral operation 
that leads to the intensity of radio-galaxies is made in presence of curved zones 
of emission. 

I express my gratitude to Miguel Onorato for useful help during the set up of 
Figure 2. 
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